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(3) Table 2 


F 
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Average flux into ground nth interval Jaeger 

(3) 


(t 

T 

P 
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Period of flux F 
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Surface temperature as periodic 

Jaeger 

Jaeger/ 

(3) 


s 

FS 
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FA 

n 

FG 
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Solar heating flux component 
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Jaeger 
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Atmospheric heating flux component 
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Ground radiation flux component 

Jaeger 
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X 
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Table 2 


A 
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Albedo 
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Wind speed 
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Rosema 
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ABSTRACT 


Thermal models from 3 sources (Watson, Outcalt and Rosema) were 
compared using similar input data and found to yield very different 
results. Watson-type models are very sensitive to albedo variations 
and to back-radiation from the sky. 

Field data with which to check the validity of models are rare and 
usually insufficient parameters were measured. Our Pisgah Crater-Lavic 
Lake data, for March 29-30, 1975, have been re-examined to indicate 
the serious discrepancy between results for thermal inertia, between 
JPL calculations (Kahle et al. 1976) and ours, when made using the 
same original data sets. 

Roof-top modelling experiments are underway to explore the practicality 
of determining thermal parameters of a known standardized material 
(Ottawa sand), when using remotely sensed data. 
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I. INTRODUCTION 

The Stanford Remote Sensing Laboratory has maintained a long and 
close relationship with thermal infrared mapping experiments. The often 
difficult acquisition of this experience has made us uniquely aware of 
both the advantages and problems of thermal mapping employing modelling 
techniques of thermal parameters with repeated diui al aircraft and/or 
satellite coverage. This same expertise with field measurement programs 
has made it painfully apparent to us that great difficulties exist in 
the determination of thermal parameters at the ground-air interface in 
an absolute sense. 

The work performed under this grant attempts to assess these 
determinations. This assessment includes a comparison of thermal mapping 
schemes, and the JPL/Stanford field effort at Pisgah Crater, California 
which indicated estimates of these "body-parameters" of the surficial 
materials may be in error as much as 100%. In light of these indications 
a study of a material with laboratory-determined thermal conductivity 
and diffusivity was begun to attempt to fix confidence limits on attempts 
to make absolute measurements of soil/rock thermal inertia. 

Adequate technology and reasonable software exists to create thermal 
inertia images from twice daily aircraft-satellite coverage, however, 
the accuracy of these determinations in terms of absolute measure remains 
to be established (Kahle, et.al., 1976). 

II. COMPARISON OF THERMAL MODELS 

The prediction of diurnal surface temperatures, given the thermal 
parameters of a material, is of considerable significance to any thermal 
study. The contemporary principles of one-dimensional heat flow on the 
earth's surface can be converted into mathematical models to resolve 
surface temperatures, and have been programmed by a number of scientists. 
The mathematical model is an extremely valuable technique, and infrared 
studies are Improved by understanding the modelling results and using them 
to help recognize variations in thermal properties and anomalous geothermal 
sites. 


In order to study and perhaps analyze the significance of thermal 
modelling the components of four competitive models that compute the 
diurnal surface temperature will be analyzed. These models represent 
advances in our understanding of heat flow over the past five years. 

All four models were written for digital computations since they are based 
on rather complicated mathematical solutions. By comparing results of 
the models the significance of various thermal parameters will be made 
clear, and the applicability of their results to delineation of surface 
materials will be apparent. 

A. Description of Thermal Models 

MODEL 1. Watson (1971a and 1971b) developed the first computer program 
of thermal modelling for interpretation of infrared images. A mathematical 
model for the diurnal surface temperature variation of the earth's surface 
was derived from the one-dimens ions 1 heat-conduction equation of Jaeger 
(1953) by applying the Laplace transform; 
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ir T 


E. V <{h-s+l . 

3=1 e 


( 1 ) 


where is the average flux into the ground in the nth interval, (y) 

is thermal inertia, T is the period of the flux, V is the surface 

P s 

temperature as a periodic step function, and <(in-s+l is a set of coefficients 

determined solely by the number of incremental steps m, in the period T^. 

The flux into the ground, F , is considered to be the result of solar 

n 

heating, FS ; atmospheric heating, FA ; and radiation from the ground. 


FG 


where 
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FS + FA - FG . 
n n n 


( 2 ) 


FS 


fS. (l-A)M(Z) Cos Z + or 


AD 


Eov . day 
n 


or 


AN 


- lov =0. 
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night (3) 


and f is cloud cover, S. is the colar constant, A is albedo, M(Z) is air mass or 
an atmospheric transmission (M(Z) ■ 1-0.2 sec Z) , Z is the zenith angle, 
o is the Stefan-Boltzmann constant, T^ is the effective nighttime sky 
temperature, T^ is the effective day time sky temperature, and 
E is the mean emmissivity. 
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The model to this point sssumes the ground Is fist, therefore, 
to improve the model a topographic correction was Introduced. 
that a change In surface slope .ill only alter significantly the incident 

solar flux. 

FS = fS. (l-A)M(Z) Cos Z'. ' 

n 

.here Cos Z* Is the zenith angle modified for local slope of dip 
angle d and strike azimuth from north <j>. 

Cos V - cos(X-dslnp) cos6(cos(.t+dco.i) + tanU-dsimtOtend) . 

and X is ih local latitude, 6 is sun's declination, t is local time from 
noon, . the diurnal angular frequency. Significant topographic variation 
.ould therefore produce, for a single site, a large variation in *« £ “« 
remperature that .ould he at a maximum near noon and at a minima near da™. 

Atmospheric effects are treated in terms of the cloud cover factors 
and transmission as a function of zenith angle. The model 
variations in convection, evaporation, and condensation in the night 

day sky temperature terms. 


night 


r.-s.-arjsss.ws sras ™ s,™ ““ 

now expresses the relation. 

OOS - coed co.Z -slndfsindcosd sin »t - co. t .indco.X - sind slnl cos .t) . 

However, because »e have t.beh ^ 

follow to be zero the results are c P J , dl r the correct 

^“d^r^r^p 1 :^ -»> - - — 1 

zenith angle relation, 

cos Z - cosX cos5 cos wt + sinAsinfi 


which is correctly formulated in the model. 


From analysis of computer runs Watson (1971a) concluded that daytime 
variations of the sky temperature have a small effect on the ground 
temperature changes because insulation is the major source of surface 
heating; while at night, however, it has a marked effect on surface 
temperature because it is the major source of heating. 

The radiation from the ground, FG is affected by a number of 

n 

surface state effects. 

FG - ZoV 4 (6) 

n s 

Equation (1) is solved subject to equations (2) through (6) by 

making an initial estimate of V and iteratively improving the solution. 

s 

The surface parameters of thermal inertia, albedo, and emissivity help 
formulate the initial guess. The properties of thermal inertia and 
albedo also strongly affect the amplitude of diurnal temperature 
variations and therefore significantly affect the accuracies of the 
model calculations. 

Watson et al. (1971b), has employed the WATEMP model in the study 
of an area near Mill Creek, Oklahoma. Model results agred with thermal 
contrasts exhibited by thermal infrared images from altitudes of 150m to 
17 km. 

MODEL 2. Subsequent to the development of the Watson thermal model a 
similar program was developed, also based on a sinusoidal input for 
solar flux, at the Stanford Remote Sensing Lab. The Stanford temperature 
model program, SURTEMP, was independently written early in 1972 by 
Andrew Green for fully interactive use on IBM 360/67 and subsequently 
adapted for use on the PDP-10 computer by Frank Honey. The program 
is intentionally similar to that of Watson in that input consists of 
cloud cover factor, albedo, emissivity, thermal inertia, latitude 
and sun's declination, orientation, and inclination of the surface, 
and the day and night sky temperatures. However, differences exist 
in the computations using thermal inertia and emissivity, and in 
establishing the initial estimate of surface temperature. The Watson 
program .teratively solves the equation for $ values until AT between the 
last two iterations is 0.1°K. SURTEMP solves the equation employing a 
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maxlmum of 20 iterations (<J> values). From the work of Jaeger (1953) this 
should yield an error in amplitude of the approximate results of about 2%. 

MODEL 3 . A digital surface-climate simulator was developed by Outcalt (1972) 
following an analog solution for the diurnal surface thermal, and energy 
transfer, regimes based upon equilibrium temperature solutions by Myrup 
(1969). Though the model was developed for climatological studies, it is 
nevertheless applicable to surface thermal infrared studies in that it 
computes the diurnal variation of temperature from the surface to a depth of 
25 cm. The basic equation of energy transfer or radiation balance, 

R - H + LE + G. (7) 


set equal to zero, is the basis of the Myrup solution. Where H is the 
vertical transfer of sensible heat to or from the air column, Le is the 
evaporation, and G is the net rate or flux of sensible heat into or out 
of the soil. By characterizing these terms by measurable parameters 


the Outcalt model takes the form: 

(l^)(Qfq) + ZoT sky 4 


4 P fc2 U 2 

EoT + 2 x 

° (ln)Z /Z )r 

L O 


C(T 2 -rZ 2 -T Q ) + L q 2 -(X w /L) f(T Q ) + 


K 

-J5 (T 

(Z g /2) U n 


■v 


0 . 


( 8 ) 


where Z 2 * k, r, P, are fixed constants: (Q*-q), T gfcy , u 2 * q 2* T 2* are 

meteorological variables; and a, E* Z q , X^, Z g , K g , are terrain variables. 

The number of boundary conditions of the Outcalt model greatly exceed the 
Watson or SURTEMP models. However, the model results contain estimates of 
climatological and soil conditions useful to the climatologist (see Table 2) . 

MODEL 4 . A most recent thermal model was developed by Rosema (1974) to calculate 
the diurnal variation of surface temperature. The program solves the heat 
transport equations numerically using parabolic differential equations. It 
initially calculates soil conductivities and heat capacities, followed by 
calculations of matrix potential and temperatures, and heat and water balance. 
Input of the Rosema model like that of Outcalt is of greater relevance to the 
climatologist or soil scientists. 


B. Comparison of Thermal Models 

In sn attempt to understand the implications of the differences in 
input parameters and mathematical solutions between these four models, data 
sets and output temperatures were compared. The SURTEMP program was already 
established on the Stanford PDP-10 system and therefore available. Watson's 
(1971a) program had been published, and after modifications to achieve 
compatability with the Stanford IBM 360/67 system by this author, it too 
became available. (Cost of the Watson thermal model (WATEMP) program on 
the Stanford IBM 360/67 far exceeded SURTEMP runs; each input set required 
approximately 0.75 minute CPU Time or roughly $6.00). Although the Outcalt 
and Rosema modeling programs themselves were not available at that time, 
test run results had been published and these results could therefore be 
compared to either program available at Stanford. 

R.J.P. Lyon (1974) over the past several years has run comparisons 
of the Friedman (1968) field data to the SURTEMP model. 

Carefully documented and measured field data relating the model 
input parameters to surface temperatures are quite rare. The only 
significant, geologically relevant, study to date is that of Friedman 
(1968), but though this study did relate surface temperatures to albedo, 
emissivity and thermal inertia, sky temperatures were not measured, and no 
calculations of thermal models appear. By assuming the soil parameters 
are correct and knowing that the cloud cover factor for those days was 
zero, the SURTEMP program was run, varying slope orientation and inclination 
and sky temperatures until model temperatures matched observed temperatures. 

By comparing At between maximum and minimum temperatures Lyon 
was able to conclude that surface dip slope changes of 15 to 50Z and 
surface strike changes of 10 to 33X, will produce a 1°C change in maximum 
or minimum target temperature. Watson (1974a) compared the differential 
change in property values which produced a 1°C change in predawn temperature 
and arrived at similar results. Other parametric variations, which will 
produce a 1°C change in target temperature are listed below; 


TABLE 2 


INPUT 

WATSON-SURTEMP 

Cloud Cover (f) 

Albedo (A) 

Emissivity (£) 

Thermal Inertia (y) 
Latitude (X) 

Suns Declination (6) 
Strike (<)>) 

Dip (d) 

Day and Night Sky 
Temperature 

<V <V 


PARAMETERS OF THE THERMAL MODELS 

OUTCALT ROSEMA 


Solar Declination (DEC) 
Radius Vector of Sun (R) 
Sky Radiant Temperature 

<W 

Dust Content (D) 

Air Temperature (T^) 

Hind Velocity (u) 

Air Humidity (q) 
Precipitable Water (w) 
Station Pressure (p) 

Soil Vol. Heat Capacity 

(GC) 

Soil Thermal Diffusivity 

(GD) 

Surface Roughness (Zo) 
Albedo (ot) 

Wet Fraction (Xw) 

Shadow Fraction (SHDRAT) 


Albedo (a) 

Emissivity (E) 

Porosity (X^) 

Vol. Qtz. Content (X ) 

m 

Vol. Clay + Feldspar Content (X ) 

m 

Vol. Organics Content (Xo) 
Latitude (<j>') 

Day (Suns Declination) (6') 
Aerodynamic Roughness (Zo) 
Windspeed (v) 

Air Temperature (Ta) 

Air Humidity (s) 

Depth of Subsoil Water Level (h) 
Temp at Subsoil Water Level (T) 
Moisture Characteristics (0 ) 
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t 




Percent Variability for 1° Change 


Parameter 

iZ2“ 

Watson 

Albedo (A) 

1.5-3 

17 

Emlssivity (£) 

25 

2 

Thermal Inertia (Y) 

12-15 

10 

Sky Temperature T or T^ 

1-2 

— 


It is apparent that a considerable discrepancy exists between the 
calculated results of Watson and Lyon, for albedo and emissivity while the 
thermal inertia values are quite similar. Watson (1974a) did not calculate 
a percent change for sky temperature variation. In an attempt to resolve 
these discrepancies, the Watson model was run using the same Mono Lake data 
(Friedman 1968) used by Lyon to form his conclusions. These comparisons 
are plotted on the following figures, input parameters and comparison of 
maximum, minimum, and temperature differences appear on each graph. 

Initial comparison involved using two different rock types. Figures 1 
and 2 compare the programs for olivine basalt icpilli (Qvb) and silty playa 
and deltaic deposits (Qal) . Inspection of the variations indicates that the 
nature of the diurnal surface temperature curves sre quite similar. Maximum 
and minimum temperatures depart from one another by only one degree for 
Qvb, and the AT for both runs are within 1°C. Comparison runs of sky temperatures, 
(Figures 3-4) , using the same input data for lacustrine carbonate (Qm) , were 
made. The shape of the curves are again similar and variations in AT are 
less than 1°C. 

Variation between the two programs appears to be unrelated to the sky 
temperature chosen. 

Comparison runs of albedo and emlssivity are compared on Figures 5-6. 

The curves are again quite similar in shape, and AT is again within 1°C. 

A number of preliminary conclusions can be made concerning these two 
models. As one would expect the programs are quite similar, and the 
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gives lower temperature estimetes for any particular hour. The shape of 
the curves are quite similar and temperature differences ere minor except 
during the late evening to predawn hours. 

A comparison of the Outcalt (1972) thermal model to the Watson 
(1971a) and SURTEMP models was made using published input and model results 
of a test run inserted both into the Stanford IBM 360 WATEMP program 
and the PDP-10 SURTEMP program. Input parameters given by Outcalt (and 
corresponding to the WATEMP-SURTEMP model) are: 
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Latitude ■ 49.3° 

Solar decl. ■ -14.9° 

Albedo ■ 0.15 
Sky radiant T. « 235°K 
Soil thermal diffus. » .0056 (cgs) 

Soil vol. heat capac. * 0.500 (cgs) 

Shadow ratio « 0.0 

From the diffusivity (a) 0 nd volumetric heat capecity (C) it is possible 
to compute conductivity (k) end thermal inertia CY) from the relation 


a 


k 

P c 


(kC) 


1/2 


0.0056 - k/0.500 


k ■ 0.0028 cal cm 
sec” 1 °C~ 1 


-2 


Y - (0.0028 x 0.500) 1/2 


y * 3.74 x 10" 2 cal cm” 2 
sec” 1 °C” 1 


An emissivity of 0.98 was assumed from the description of the site-"Needle- 
Ice Event", and the slope orientation and inclination are assumed to be zero 
from the shadow ratio. Additional input for the Outcalt model included values 
for air column and soil conditions. 

Results of the comparison run are illustrated on Figures 7 and 8. The 
shape of both curves are similar, and maximum temperature varies by only 1°C, 
minimum by 4 C, and AT by 3°C. Interestingly, again the largest discrepancy 
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between the two models exists In the predawn hours. From the comparison 
It appears we can tentatively conclude that both programs yield similar 
temperature estimates (^4°C) . The additional complexity of Outcalt's program 
however does not necessarily make its results more accurate. Uncertainty 
of parameters must still exist and added parameters may increase total 
uncertainties. To ensure the usefulness of a model the additional parameters 
that are of value only to the climatologist should be disregarded by the 
geologist or geophysicist. This would make the job of collecting field 
data simpler and the model less cumbersome to use. 

Two comparisons of the ROSEMA/WATEMP / SURTEMP models were made using 
the input parameters published by Rosema (1974) for coarse sand and Basin 
clay at sites in Holland (52°N), that he does not describe in detail. Input 
values given by Rosema applicable to the Watson program are: 


Coarse Sand Basin Clay 


Albedo 

0.30 

0.15 

Emissivity 

0.92 

0.94 

Vol. Qtz. Content 1 

„ 0.60 

0.00 

Vol. Clay + Feldspar) 

m 0.00 

0.46 

Vol. organics (X ) 

0.00 

0.00 

o 

Latitude 

52.0° 

52.0 

Day 

April 12 

April 


Values of sky temperature and thermal inertia had to be estimated. 

The selection of sky temperature for the location and season of the test* 
at a ground temperature of 280°K was fixed somewhat arbitrarily at 260°K 
night and ^ 7 0°K day.*** Selection of thermal inertia was somewhat less 
arbitrary. Employing the De Vries (1933) equation for soil heat capacity: 

C - 0.46X a + 0.60 X q + X y cal cm -3 °c" ' (9) 

the coarse sand with a moisture content of 0.20* would have a heat capacity 
of: 

C - 0.46(0.60) + 0.60 (0.0) + (0.20) - 0.476 (10) 


*** 


Humidity at the site should produce close night and day temperatures 
(Lyon* personal com mun ication 1973). 


with a moisture content of 0.20; the Basin clay would have a heat capacity of: 


C - 0.46 (0.46) + 0.60 (0.0) + (0.30) - 0.512 (11) 


Estimates of conductivities given by Sellers (1965) and De Vries (1963) 
at moisture contents of 0.20 are 4.2 x 10 ^ and 5.3 x 10 cal cm sec C 

for coarse sand; and 3.5 x 10 and 2.8 x 10 cal cm sec C for Basin 
clay. Taking the mean of these values with the calculation of estimated 
volumetric heat capacities calculated , values of thermal inertia (y) can 
be estimated: 


Y 

Coarse sand 


p (kC) 


//. a 


1/2 


.“3 


a i n\ 


1/2 


( 12 ) 


- 0.048. 

_-i 1/2 

Basin clay * (3.2 x 10 x 0.51) 

- 0.040. (13) 


Results from the comparisons (figure 9) indicate a large discrepancy of 
surface temperature. The comparison, nevertheless, was a useful exercise 
in showing the difficulty of obtaining thermal parameters for any given 
test. The large uncertainty in our estimation of appropriate sky temperature 
and thermal inertia illustrate their significant effects on modeling of the 
diurnal temperature variation. 

Certain added ambiguities exist in the Watson-Rosema Comparison: 

Within the description of the parameters of the sample materials, by Rosema 
(1974) , it is difficult to imagine an absolute zero quartz content in the 
Basin clay. Also somewhat surprising are the results given by WATEMP. It 
is unusual that surface temperatures in early April at a north latitude of 
52° could reach 50°C. 


C. Discussion 

The value of thermal models seems apparent from the previous discussion. 
The thermal models of Watson and SURTEMP eppear quite similar, and hold 
an advantage over those of Outcalt and Rosema in their simplicity of input 
parameters. When used to differentiate rock and soil types as on infrared 
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^sery, the model Indicates that thermal inertia is a good discriminator 
of materials, and the effects of thermal inertia on surface temperatures are 
maximized near dawn. Thus dawn represents the optimum time to observe 
thermal contrasts due to thermal property differences. 

HI* PISGAH-LAVIC LAKE STUDY (March 29-30, 1975) 

A. Site Description 

The Pisgah Crater-Lavic Lake area is located in the Mohave Desert, 

San Bernadino County, California. It is approximately 38 miles east-southeast 
of the town of Barstow. Pisgah Crater and its associated flow is nearly 14 
miles long and 4 miles across. At the southeastward end of the flow is Lavic 

Dry Lake approximately 2 miles in diameter. Other details of the area appear 
in Kahle, et.al., 1976. 

The Pisgah Crater flow is composed of numerous thin olivine basalt flows. 
The top of this flow sequence is made up of both pahoehoe (ropy) and aa (jagged 
clinkery) lava. The aa lava is dark gray, or more usually, black in color and 
has an extremely rough surface. The pahoehoe lava is a "medium" dark gray 
in color and has a smoother undulatory surface. Both types of surface 
can accumulate windblown sand or silt. Pisga., Crater is ,aade of cinders 

(that range in color from black to grayish red), and small bombs. The crater 
floor is pahoehoe basalt. 

Lavic Dry Lake is a very flat playa composed of a very pale, yellowish- 
brown silty clay. Neal (1965) determined the playa is made up of 79Z clay 
(essentially montmorillonite and illite) 20.7Z granular components, and 0.2Z 
accessories with a bulk density of 1.67. The playa had an approximate maximum 
relief of 2 1/2 feet. Numerous areas on the playa are mantled by pebbles or 
cobbles of basalt or alluvial fan material. These mantling materials are 
somewhat anomalous as they are found some distance from their obvious source. 
Gawarecki (1964) reports they are believed to be ice rafted to their current 
position when sufficient lake levels existed. 

Alluvial fan material is found along the borders of Lavic Lake. This 
material is a heterogeneous conglomerate derived from present topographic 
highs. It is a yellowish gray color with a sparse vegetation cover. 
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Measurement sites at the crater, designated, "Crater outeredpe " ‘ and 
"C rater station" consisted of cinders, clinkers and bombs; site "Crater outside 
Basalt" consisted of aa and pahoehoe basalt east of the cone viewed from a 
distance. Site "Crater center" was measured from a distance and consisted 
of solidified pahoehoe lava; Site "Crater Red Cinder" consisted of an area 
of predominantly Red cinder on the inner lip of the crater; and site "Crater 
8oad"consisted of compacted cinders and dirt. 

Lavic Lake sites 1-1, i_2 i_-j i_a 9 ■> o •s . , . 

■*■*** 2-2, 2-3, consist of the pale 

yellowish-brown silty clay. Site 2-1 was on alluvium (Figure Map 1). 

B. Albedo Measurements and Analysis 

Measurements of albedo were made at the test sites on Lavic Lake and at 
the crater employing an ISCO spectral radiometer and an EXOTECH radiometer 
by Stanford Remote Sensing Lab personnel. These measurements were necessary 
to better understand the site locations and to provide the necessary input data 
for our thermal model. The procedure involved making a number (5-25) of readings 
at each site attempting to include the greatest diversity in surficial 
conditions, with concurrent readings of white reflectance standards of BaSO 
and Fiberfrax. ISCO readings (see Appendix 2) were made from 0.400 to 1.55* 
micrometers with a bandpass interval of .025 micrometers in the visible 
and a bandpass interval of .050 micrometers in the near-IR. Exotech readings 
were collected in detail at the Lavic Lake sites. The instrument sensed in 
the same channels as the LANDSAT MSS sensor (0.5-.6, .6-. 7, .7-. 8, .8-1.1 
micrometers). Twenty-five readings at Stations 2-1 (alluvium), 1-1 and Base 
camp (playa) were made, and five readings at Stations 1-2, 1-3, 2-2A, 2-2B, 

2-2A (playa) were made (Figure Map 1). Appendix 1 catalogues these measurements 
and yields values of bandpass, reflectance, and'pseudo-ClE^cclor coordinates 
for each channel and calculates an albedo for each site. 

The albedo values from the Exotech measurement sites are compiled in 
Table 3 along with means and standard deviations for each location. The 75 
playa measurements yielded a mean of 0.375 and a standard deviation of 0.027 
(7% variation) while the 25 alluvium sites (2-3) yielded a mean of .295 and a 
standard deviation ot .018 (6Z variations). These statistical checks indicate 
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TABLE 3 


Albedo Determinations i 

PLAYA 


Site* 

1-1 
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V 

A total 
s total = 
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.009 

Sitei 

1-2 

1-3 

2-2A 

2-2B 
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,4o4 

.365 

.385 

.461 


.365 

.369 
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.386 

.469 


.373 

.375 

.348 

.388 

.457 
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.372 
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.375 

_a2®4 . 
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Table 3 continued 


ALLUVIUM 


Site* 

2-1 

2-1N 

2-1E 

2-lS 

2-lW 


.293 

.299 

.288 

.305 

.298 


.28? 

.273 

.358 

.290 

.298 


.293 

.282 

.278 

.300 

.272 


.290 

.291 

.311 

.303 

.271 


. 286 

.288 

.*22§_ 

tin 

.. £12 

X= 

.290 

• 

r\> 

00 

^3 

.307 

.302 

.290 

s= 

.003 

.010 

.031 

.008 

.018 


X total = * 295 
s total * ’° 18 


n = 25 


Emissivity for the playa was assumed, at each site, to 
be .985. 
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— " th ‘ ^ tautochronea ,or depth-temperature plots. 

«.la of th. tautochrona. permit poaitioning to within ,5°C and .5 c 

HI Ch ' ref0re ‘ hOUld b ' —at. rapreaantationa. Append!, 

lineT T 8 “ te<l (aV ‘ raS ‘ ° f 5 r “ dln8,) ‘ Urf «* temperatures (.olid 
line) compared with runs of the SURTEMP and WATEMP programs 

Th. najor concarn in the determination of input paranatal, ant.il. 

aurflceT 10 " ‘ V * 1 “‘ ! (Y) 60 “«*»taly natch tha meMured 

ourfaca temperature.. Varloua nathod. wara -ployed to arri.a at a 

hat can „a conaidarad rellahla. Input parameter. that arc avidant for 

the dry lake sites, are; 

Latitude: 34.66° 

Dip : 0.0 

Strike : 0.0 

Sun's Declination: 03°02.9'. 

The work of Neal (1965) provided pertinent data on tha bulk denaity 

rrrrr the Lavic uke piayai “ u — -*» -» 

ng tha niaaion a„d provided d . e . yleldlng uelght , ^ a( „„ ^ 

ract on .20Z). With thia data it is possible to apply the DeVries (1963) 
oruula for volumetric heat capacity (C) of the play, notarial. 


c - 0.46 x + 0.60 X + x , 
m o w 


(14) 


where X, - vol. tract, ninar.la X g . vol. tract, organic, and X - vol. tract 
water DeVries ala. detained a fornula for th. danping dap^h of . .on ' 
elated to thermal dif fusivity ( a) , 


(2 a / w ) 1/2 . D 


(15) 


where D 


angular period of the diurnal wave 


Damping depth, and w 

t™ 76 ' 10 " 5 r ‘ d ‘ eC ' 1) • ^ d “ Pl "» d * pth “ depth wh.r. th, aurf.c. 

inspH!i Ure ^“ idered at “ 8lnusoid *l wave, is reduced by 1/e (.37). From 
paction of th. tautochrona. at th. play. aita. th. dating depth occur, 
at xu cm. 

.« . rzrrr " c,n ,oiva f ° r tbanui diftu * i ’ u >' — •«!». 

X cm sec . if we consider this value and that 
or volunatric heat capacity (0.568) to b« valid than th. th.m .1 in.rtl. 
th. play, can b. c.lculatad af tar fir.t datamining th.m .1 conductivity (K 
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Now, , 

a - k/C 3.64 x 10 - k/.568 

k - 2.07 x 10 -3 cal cm -1 sec -1 °C _1 (16) 

y - (kC) 1/2 - (2.07 x 10 -3 x .568) 1/2 

■ .034 cal cm 2 °C ^sec . (17) 

Additional information needed for the SURTEMP thermal model included 
a cloud cover factor assumed initially to be zero, a hemispherical net sky 
radiation term (sky temperature) for the day and night, and the emissivity 
of the surface approximated from measurements of similar surfaces at 0.985. 

From previous work (see Section II) of Lyon (1975) we note that an inaccuracy 
in the emissivity of 252 will produce only a 1°C change in surface temperature, 
employing the SURTEMP program. However, as we have noted, the thermal model 
is quite sensitive to sky temperature. Unfortunately a net radiometer, 
necessary for these measurements, was not available and thus sky temperature 
could be only approximated from the known meteorological conditions. These 
calculations and approximations along with the albedo measurements were 
performed for each of the playa sites and dlosestfits to the observed data 
were plotted (Appendix 4). As can be seen, fairly close fits can be achieved 
with day and night sky temperatures of 240°, 230°K respectively and employing 
the calculate thermal inertia (Y) of .034 cal cm °C sec . 

Discussion of these results with the JPL investigators, whom we were 

assisting, indicated that they believed a more accurate thermal inertia value 

for the playa fell between .01-. 03 cal cm 2 °C \ The thermal model developed 

by Kahle (1976) compensates for wind at the site and they felt this 

meteorological variable (not directly included in the SURTEMP model) invalidated 

our results. To further explore this discussion additional SURTEMP runs were 

made analyzing station 1-1, (Appendix 5). It was determined by using the cloud 

cover factor to simulate a moderate wind speed (i.e. loss of input heating) 

and by altering the sky temperature the results could be made to approximately 

—2 o —1 —1 

conform to a thermal inertia of .02 cal cm C sec . 

Unfortunately, we have not convincingly proven, even in the simple case 
of the playa material , what is the true value for the thermal inertia! Though 
admittedly the SURTEMP program does not include a compensating factor for 
strong windspeeds, it must also be pointed out that the JPL results were based 
on an assumption that the playa moisture content was zero! It can be shown 
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that this assumption can reduce a final thermal Inertia value by more than 
30%. Scientific integrity would therefore force us to conclude that an 
exact measure of thermal inertia was not determined by either group. If it 
is to be accurately estimated, precise measures of net sky radiation, emissivity, 
soil moisture, and wind speed must be made and included in thermal modeling. 

It is abundantly clear that solution of thermal models using real data 
is not a simple problem, and much remains to be done. The problem is 
especially complex if the measurements are made from a remote vehicle. 

IV. SANDBOX STUDY 

A. Objectives 

In an attempt to determine the accuracy of thermal inertia determinations 
by remote sensing techniques a study was begun at the Remote Sensing Lab. 

Though many workers are now reporting thermal Inertia values using various 
computer model deviations of the heat conduction equation, it is difficult 
to fix confidence limits to these determinations. Therefore, after a considerable 
search, a standard material was found for which both thermal conductivity and 
diffusivity had been determined by direct laboratory techniques at varying 
moisture contents. A heat transfer study by Moench (1969) had determined these 
parameters for 20/30 mesh Ottawa sand (quartz, 0.5 to 0.8 mm size). 

The sand could thus be placed under simulated natural conditions and attempts 
made to correctly determine its thermal properties. If successful further 
Impetus would be gained for heat capacity mapping. If, however, a thorough 
experiment proved completely unsuccessful, we would be forced to examine the 
total validity of thermal Inertia determinations and their discrimination 
by remote sensing techniques. 

We have only begun preliminary experimentations, but these results 
and tentative assessments will be reported. 

B. Measurements 

In order to simulate natural soil conditions the sand was placed within 
a square box, constructed of Douglas Fir and sealed with water-safe glue, 

screws. The thickness of the walls is uniformly 6.3 cm, and it is 
30.5 cm deep and 45.7 cm long and wide. The sand and thermistor probes 
placed within the sand were allowed to equilibrate for over one week. 
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For thi8 preliminary work the sane infomation was collected as 
was collected during the Plsgah-Lavic Lake mission. This procedure was 
adopted In order to assess the importance of accurate, in situ, net sky 
radiation and emissivity measurements. All other parameters were accurately 
measured, including albedo (.5-l.lUm) which was determined to be 0.54. me 
thermal measurements are compiled in Table 4 and the tautochrone of the data 
is depicted in figure 12. The emissivity of the sand was considered to be 
about 0.95, and the sky temperature (day and night) were both estimated 
at 263°K. Partly cloudly conditions existed sporadically during the day 
and evening and seemed to indicate a slight cloud cover factor. 

C. Results 

The water content of the sand was kept as close to aero as possible, 
yielding, from Moench's work, an expected thermal inertia of 0.0166 cal 
cm C sec . Figure 13 represents the best fit for the observed data 
with these input parameters, me 5°C discrepancy between the measured and 
modeled final temperatures, when the thermal inertia is known, succinctly 
states the need for more accurate measurements of net sky radiation, emissivity, 
and cloud cover-vindspeed I 

During this NASA grant a "look-up table" program for thermal inertia, 
based uponSURTEMP , was developed, me program allows for the construction 
of a 3-dimensional matrix, developed from the standard (or fixed) input data 
sets and with variable albedo and day and night temperature as the sand 
ordinates. Thus for any particular time of day and or night, a similar table, 
as seen in figure 14, can be generated. In this instance, the program yields 
a thermal Inertia value of 0.0195 cal cm 2 °C l sec"’ 1 . me input is identical 
to that of figure 13. Comparison of these techniques allows us therefore to 
conclude the inaccuracies in our input parameters has produced a discrepancy 
of .003 in the thermal inertia. 

These efforts are continuing and will be expanded with renewed 
funding for the HCMM program. 

V. CONCLUSIONS 

A. Little agreement exists today between results from several thermal models 
using the same input data sets. 

B. Each model has a varying degree of sensitivity to any specified parameter. 
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mostly to sky temperature. 

A net radiometer appears to be a vital part of any ground installation, 
and the net-radiation balances as a function of time of day and 
weather, are necessary. 

Fitting remotely-sensed data to spot locations around weather station 

— teS > 011 flaC 8round ’ may be possible, but extrapolation to outlying 
areas is unlikely. 

Relative determinations of thermal properties may be possible using 
remote sensing, but absolute value determinations appear most unlikely. 
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APPENDIX 2 


ISCO Bi-DIRECTIONAL REFLECTANCE MEASUREMENTS 


^ •, f> . w-v. „ J 
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DATE: March 28, 1975 

TIME: 1205 

LOCATION: Lavic Lake 

• 

STANDARD : BaSC> 4 

OBSERVER: Gary Ballev/ Stuart Marsh 


TARGETS: #1 - 1-1B 

if 2 - Base Camp 

if 3 - Alluvium 


Time: 1131-1210 

1416-1431 

1433-1459 1555- 

• 

Wavelength 


7*0 OTT1 


■Ill 


1.25 

1.75 

0 

3.5 

10.0 

1 

11.0 

38 

6 

17.5 

55 

9 

18.9 

57 

10 

19.4 

56 

10 

19.8 

52 

10 

19.5 

47 

9 

31.0 

66 

15 

33.0 

66 

16 
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DATE: 

March 29, 1975 LOCATION: 

Pisgah Crater 




ISCO HEAD: 

Probe OBSERVER: 

G. Ballew 



TARGETS: #1 - Cinders Near W. 

Cone 





in - Red 

Cinders, 20 

' West of 

n 



TIME: 

1036-1054 

1555-110 

118-1133 1347-1400 

1310 

Wavelength 


BaSo4 

Fiberfrax 



Bi 

(ym) 



SA = 51° 

SA = 

53° SA = 54.3 

O 

.400 


1.80 

4.90 

0.24 

0.30 

1.55 

.425 


8.0 

10.5 

0.67 

0.80 

8.10 

.450 


29.0 

34.0 

2.45 

2.65 

27.5 

.475 


42.0 

50.0 

3.50 

3.80 

38.0 

.500 


43.0 

50.0 

3.60 

4.10 

40.0 

.525 


41.0 

48.0 

3.90 

4.30 

39.0 

.550 


40.0 

45.0 

3.95 

4.25 

31.5 

.575 


34.0 

38.0 

3.80 

4.20 

33.0 

.600 


51.0 

56.0 

5.95 

7.00 

47.0 

.625 


52.0 

56.0 

6.10 

7.50 

48.0 

.650 


51.0 

55.0 

6.10 

7.85 

47.5 

.675 


45.0 

49.5 

5.65 

7.40 

42.5 

.700 


42.0 

45.5 

5.40 

7.25 

39.0 

.725 


37.5 

40.5 

4.85 

6.70 

35.0 

.750 


31.5 

34.0 

4.10 

5.90 

29.0 

.75 


51.0 

55.0 

6.30 

8.20 

43.0 

.80 


50.0 

52.0 

6.10 

8.50 

37.0 

.85 


41.0 

46.0 

5.30 

2.25 

21.0 

.90 


32.0 

35.0 

4.00 

5.70 

21.5 

.95 


24.5 

26.5 

2.95 

4.30 

23.5 

1.00 


26.5 

29.5 

3.10 

4.90 

22.0 

1.05 


24.5 

27.0 

2.90 

4.75 

20.0 

1.10 


22.5 

23.5 

2.65 

4.50 

17.3 

1.15 


9.0 

10.5 

1.20 

2.00 

7.2 

1.20 


9.0 

9.5 

1.15 

2.15 

7.6 

1.25 


16.5 

18.5 

2.15 

4.20 

14.0 

1.30 


18.5 

21.5 

2.45 

4.90 

15.5 

1.35 


10.5 

12.0 

1.25 

2.35 

8.0 

1.40 


2.25 

2.50 

0.29 

0.30 

1.70 

1.45 


1.50 

1.80 

0.21 

0.40 

1.15 

1.50 


4.10 

4.70 

0.55 

1.20 

3.4 

1.55 


6.00 

7.10 

0.83 

1.75 

5.0 
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MEASURED SURFACE AND PROBE TEMPERATURES 
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42 . 0 

43.0 

43.0 

42.0 
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38.0 

38.0 

67. 

FR 20 
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30. 0 

55.0 
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38.0 

69. 

PS 20 
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32 .0 
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69 . 
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71. 

CS25 
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.742 

• 74 7 

.746 

.752 

.753 

0.5 

1.0 

0.6 

1 .0 

1 .0 

72. 

CE2 '» 

0634L 

.732 

.744 

.732 

.735 

.740 

0.4 

0.6 

0 . 2 

0.2 

0.4 

73. 

f*E<5 
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.311 

.344 
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.243 

10.5 
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8.6 

8 .R 
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11.0 11.0 
13.5 10.5 


9.2 


i . ; .. *<•>. it. 5 io.9 

, 620 16,5 14.0 11.5 


2-129 2030 f.C 13.5 
2-125 2300 3.7 f. c" 

2-130 0120 1.0 6.6 

2-130 0345 -l.o <.2 
2-120 OoOO -2. a ;. £ 
2-130 0501 c.2 *.c 

2-120 1 1C 3 24.6 12. f 
2-1:0 134) ?3.4 ; 0 .1 
2-120 1500 
C IA\IC L^K 
2-225 0502 


1222 2 


13.5 12.0 
11.2 11. 9 
11.0 12.0 

9.7 !<.. 0 

7.7 H.9 
7.9 12.3 
3.6 12.0 

11. 1 ll.O 


15.6 


24.4 

22. * 

13*1 ll.O 


2 2-2 
-1.3 

2.0 

t-nObi 
8.1 11.5 

18.1 

-1 . 2 

l.E 

7.4 11.5 

13.0 

22*7 

t.e 

10.4 H. j 

13.1 

23.5 

25.5 

13.3 
16. C 

8.6 11.3 

10.1 11.0 

15.4 

15.0 

22.0 

16. 5 

11.6 10. 4 


x 1*Z * ni: ’ ll.r 13.0 il.o 15.5 

2-*t1 2iC5 4.2 1C. £ 13. u il.5 * 


12.9 



• 

V 



V 



183. 

2-2 25 

2251 1.3 

7.E 

12.t 

11. 9 


12.7 

m. 

2-2 *f> 

0155 -1.7 

A. 5 

10.3 

xi. 7 


14 . . 7 

165. 

2-2 ;C 

0-35 -2.0 

1.6 

6.0 

n.6 


ti.O 

laf . 

2-220 

0120 -t.O 

C.C 

7.5 

X 1 . 1 


12.5 

197. 

:-2:-o 

0518 ID. 6 

*.2 

7.2 

It. 9 


lj.3 

IB 1 *. 

2-230 

1-1 5 2 2’. 9 

12.3 

7.o 

11.8 


13.2 

109. 

2-220 

1330 32.5 

20.1 

10. u 

10. 8 


U.5 

190. 

2-220 

1635 32.5 

22.1 

11.1 

10. 5 


13.1 

151 . 

C LAV IC 

LAKE 2-3 


PKUdE 



!91 . 

2-32° 

0506 -2.0 

:-.! 

e .2 

It. 9 



1^3. 

2-225 

0 7 2*» -2.0 

1.5 

7.0 

10.6 



166. 

2-2 25 

1012 K.O 

6.C 

C • 9 

it. a 



155. 

2-2 25 

1225 26.0 

12.1 

8.*. 

11.0 



15o. 

2-225 

1623 25.6 

16.6 

9.9 

10. 7 



197. 

2-225 

1613 21.2 

17.5 

11.2 

9.1 



158. 

2-225 

1818 12.2 

16.1 

12.x 

10.2 



150. 

2-225 

2 IGfc 2.0 

10.5 

12.0 

11.0 



200 . 

2-22S 

2253 -1.0 

E.C 

11.8 

11.2 



201. 

2-220 

0116 -6.3 

5.0 

10.2 

11.0 



202 . 

2-220 

0336 -5.5 

6.5 

6.0 

il.O 



203 . 

2-320 

0621 -B.0 

1.0 

7.8 

it.O 



2‘J*». 

2-220 

0517 15.6 

6 . ( 

7.5 

11.6 



205. 

2-2 20 

1023 31.6 

11. C 

8.0 

li.6 



206 . 

2-2 20 

1321 22.2 

1 6. C 

10.1 

10.8 



207. 

2-220 

163t 25.0 

2C.5 

11.0 

10. 1 



208 . 

C LA V !C LAKE 1-6 


•PRUBc 



209 . 

1-625 

0623 -6.0 

0.0 

6.0 

11.6 



21u . 

1-625 

C 7 37 2.0 

2.0 

11.9 

12.0 



211 . 

l-*25 

1025 JS. 5 

E.5 

8.x 

12. 1 



212. 

l-ti6 

1253 26.0 

15.1 

6.5 

11. 0 



213. 

1-625 

16*5 25.8 

19.0 

12.0 

11.0 



215. 

1-626 

1625 25.0 

21.5 

16.5 

12.0 



215. 

1-625 

ie23 16.6 

17.'. 

16. a 

11.1 



216. 

1 _* 2 c 

2113 

10.7 

13.5 

11.6 



217. 

1-626 

23^6 3.0 

E.C 

12.0 

12.0 



216. 

1-6 20 

0125 -0.1 

5.0 

10. * 

11.6 



219. 

1-6 20 

0351 -2.3 

2.7 

9.2 

11.8 



220. 

1-620 

0627 -6.0 

1.2 

8.0 

11.6 



221 . 

1-630 

0520 13.6 

6.C 

7.6 

12.0 



222. 

1-6 20 

1113 27.6 

16.1 

9.0 

11. 8 



223. 

1-620 

1366 23.0 

21.5 

12.0 

11.0 



226 . 

1-620 

1666 36.2 

22.1 

13.7 

ti.O 



225. 

C LAVIC LAKE 1-5 





226. 

1-526 

(.6 20 0.0 

2 • f 

6.0 

tO.O 

11.5 


227. 

1-526 

0760 0.6 

2.1 

5.0 

8.5 

10.1 


220. 

1-525 

1020 11.0 

6 . E 

5.2 

7.6 



226. 

1-525 

1250 26.5 

15. E 

6.5 

6. 2 

9.1 


2 3*1. 

1-525 

1660 26.0 

20.6 

13.0 

6.6 

9.5 


231. 

1-526 

1620 30.0 

22.6 

15. S 

11.0 

9.8 


232. 

1-526 

1322 22.2 

21.2 

18.2 

12.7 

10.2 


233. 

1-526 

2112 3.0 

12. E 

15.0 

la. 5 

12.5 


236. 

1-5 2 C 

2303 6.5 

6.5 

13. ■* 

13.9 

12.2 


235. 

1-520 

0121 2.2 

6.5 

10.5 

12. 2 

12.0 


256 . 

1-520 

0 3 5U 1.0 

«.C 

8.5 

11. 7 

11.8 


237. 

1-530 

0636 -1.5 

2.0 

7.7 

11.2 

12.6 


233. 

1-5 30 

0523 8.3 

6.C 

6.6 

9.8 

11.6 


239. 

1-5 20 

1111 23.0 

16.5 

8.6 

9.6 

11.3 


260 . 

1-520 

1365 35.0 

22. C 

16.0 

10.0 

10.0 


261. 

1-5 20 

1665 28.0 

25.8 

11.1 

11.0 

10.1 


262. 

C LAVIC LAKE 1-6 


PKCftc 



263. 

1-425 

0625 0.0 

2.1 

7.5 

la.O 











2 44 . 

1-725 

07,7 

-1 .0 

2-o. 

1-625 

i;::t 

12.6 

<48. 

1-625 

1222 

21.2 

24-». 

1-6 2 9 

14.3 

22.2 

: 4 . 

1-625 

16 30 

2 2.0 

245. 

1-t 2 C 

l ‘25 

16.1 

2 50. 

1-6 25 

21 14 

11.5 

<51. 

1-6 2 5 

2 30 5 

7 . C 

28? . 

1-6 70 

0127 

4.8 

2'3. 

1-62C 

03 c 2 

1.0 

2 54. 

1-6 20 

C6 < p 

0 .0 

2 i 5 . 

1-6 20 

0623 

10.0 

2 r 4. 

1-6 '• C 

1112 

21.1 

257. 

1-6 30 

1348 

25.0 

253. 

1-6 3C 

14 50 

21.0 

25°. 

C UUC L15R 

1-1 

26U . 

1-125 

04*. 2 L 

. 3e0 

261. 

1-125 

06 101 

.372 

<62 . 

1-125 

0720 


267 • 

1-125 

05656 

.471 

2 >4. 

1-125 

1 2 0 7 p 

• 7 7 5 

26' . 

1-125 

I 415 P 

• 572 

2c6. 

1-125 

1557V 

.722 

2^7. 

1-125 

1805V 

.42 l 

2cS. 

1-125 

2C53L 

.786 

265. 

1-125 

22361 

.722 

27u . 

1-120 

03.001 

.67 5 

271. 

1-120 

03151 

. 6 7 6 

2 72 . 

1-123 

06 1C L 

.'52 

273 . 

1-120 

13166 

.275 

274. 

1-1:0 

1425V 

.250 

2 75 . 

C Lf\!( 

: lip e 

1-2 

2 76. 

1-225 

50 1, 

. 262 

277. 

1-225 

0S15L 

.3 72 

278. 

1-225 

0 7 70 


2 79. 

1-225 

vnov 

.475 

2 60 • 

1-225 

17J.5P, 

.774 

281. 

1-225 

1 2?r,p, 

> £2 6 

<82 . 

1-225 

14 2.JP 

.905 

2 8? . 

l-22 c 

16028, 

.635 

2 p**. 

1-225 

18377* , 

.421 

285 . 

1-225 

2 3 ” 5 L . 

.741 

286. 

1-225 

2745L, 

.734 

287. 

1-2 ?C 

0 1 0. 5 L . 

, 7C7 

287. 

1-2 2C 

02201. 

.6 55 

2 - c . 

1-2 20 

OS 1 1L. 

.611 

250. 

1-270 

1221.6 . 

.2 72 

241 . 

1-220 

14?. 7 P . 

.247 

252. 

C UUC 

LskP 

PC It 

253 . 

PC25 

0455L. 

27C 

254. 

PC29 

Os 1 31 . 

354 

29S . 

B 0 2 5 

0(35 


254. 

PC2 C 

1034 * . 

417 

2 7 7. 

BC25 

12177'. 

345 

258 . 

RC25 

1422 P. 

P 70 

<59. 

RC25 

ibTpV. 

682 , 

300. 

B C25 

180CV. 

456 

201 . 

C LA VIC 

IC<{ 

1-3 

202. 

1-225 1 

0458L. 

290 

303. 

1-129 0725 


304. 

i-2<« 

1306P. 

4 22 . 


C.C 4.6 liJ.U 


0.0 



5.4 0.9 13.1 


5.0 



16. C 11. 0 1 < • 6 


8 . j 



15.5. 13,3 t ^ . 4 


ll.l 



7C.E 13.6 12.9 


13 . 3 
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12 . b 



14. C 13.2 i< . 5 


4.4 



1C.1 13. < <1.5 


1.7 



6.2 13.o lo . 0 


- 1.1 



6.4 10.2 It, U 
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2.5 o , 0 |a. L 


-6.7 
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PKi-6 SPECIAL (Spectral Passband, 8- 14^) 


LO RANGE 


■Sj:5 

MED RANGE HI RANGE 


Meter 

Recorder 

Meter 

Recorder 

Meter 

Recorder 

(°C) 

(MV D-C) 

CC) 

(MV D-C) 

<°C) 

(MV D-C) 

- 40.0 

0*0 

0.0 

0.0 

25.0 

0.0 

- 39.0 

14.0 

1.0 

28.6 

26.0 

21.3 

- 39.0 

28.3 

2.0 

37.6 

27.0 

42.7 

- 37.0 

42.6 

3.0 

86.9 

28.0 

64.3 

- 36.0 

57.3 

4.0 

116.5 

29.0 

86.2 

- 33.0 

72.4 

5.0 

146.4 

30.0 

108.2 

- 34.0 

37.6 

6.0 

176.6 

31.0 

130.4 

- 33.0 

103.0 

7.0 

207.2 

32.0 

152.8 

- 32.0 

119.7 

8.0 

238.1 

33.0 

175.4 

- 31.0 

134.3 

9.0 

269.3 

34.0 

198.2 

- 30.0 

130.6 

10.0 

300.8 

35.0 

221.2 

- 29.0 

167.0 

11.0 

332.7 

36.0 

244.3 

- 28.0 

199.6 

12.0 

364.8 

37.0 

267.7 

- 27.0 

200.4 

13.0 

397.3 

38.0 

291.2 

- 26.0 

217.3 

14.0 

430.2 

39.0 

315.0 

- 23.0 

234.8 

13.0 

463 . 3 - 

40.0 

338.9 

- 24.0 

252.3 

16.0 

496.8 

41.0 

363.1 

- 23.0 

270.1 

17.0 

530.6 

42.0 

387.4 

- 22.0 

288.1 

18.0 

564.7 

43.0 

411.9 

- 21.0 

306.4 

19.0 

399.2 

4440 - 

436.6 

- 20.0 

324.9 

20.0 

634.0 

45.0 

461.5 

- 19.0 

343.7 

21.0 

669.1 

46.0 

486.6 

- 18.0 

362.7 

22.0 

704.6 

47.0 

511.9 

- 17.0 

382.0 

23.0 

740.3 

48.0 

537.4 

- 16.0 

401.5 

24.0 

776.4 

49 . 0 * 

563.1 

- 13.0 

421.3 

23.0 

812.9 

50.0 

588.9 

- 14.0 

441.3 

26.0 

849.6 

51.0 

615.0 

- 13.0 

461.6 

27.0 

886.7 

52.0 

641.2 

- 12.0 

482.2 

28.0 

924.2 

53.0 

667.7 

- 11.0 

503.0 

29.0 

961.9 

54.0 

694.3 

- 10.0 

324.0 

30.0 

1000.0 

35.0 

721.2 

- 9.0 

543.3 


' 

56.0 

748.2 

- 8.0 

566.9 



57.0 

773.4 

- 7.0 

588.7 



56.0 

802.8 

- 6.0 

610.8 



59.0 

830.4 

- 5.0 

633.1 



60.0 

838.2 

- 4.0 

653.7 



61.0 

886.21 

- 3.0 

678.6 



62.0 

914.3 

- 2.0 

701.7 



63.0 

942.7 

- 1.0 

725.1 



64.0 

971.3 

• 0.0 

748.8 



65.0 

1000.0 

1.0 

772.7 





2.0 

796.9 





3.0 

4.0 

3.0 

6.0 
7*0 
8.0 

821.3 
846.1 
871.0 

896.3 
• 21*1 

... - 

<.'» » A \ 

„ 



- 


~ 

w m » m m 

947.6 





9.0 

973.7 

- 

. 

- 

— 

10.0 

1000.0 






* 


High-level output. Low level output is l/20th of listed value 


Figure 2-2 RECORDER-OUTPUT CALIBRATION* 
(Standard Spectral Passband, 8 to 14p) 


LOW RANGE 


Meter 

Recorder 

(°C) 

(mv d-c) 

-20 

0.0 

-19 

24.0 

-18 

48.4 

-17 

73.1 

-16 

98.2 

-15 

123.5 

-14 

149.1 

-13 

175.0 

-12 

201.4 

-11 

228.0 

-10 

255.0 

- 9 

282.3 

- 8 

319.9 

- 7 

337.9 

- 6 

366.1 

- 5 

394.7 

- 4 

423.5 

- 3 

452.9 

- 2 

482.4 

- 1 

512.4 

0 

542.6 

+ 1 

573.1 

2 

604.1 

3 

643.7 

4 

665.3 

5 

695.2 

6 

727.1 

7 

758.7 

8 

790.5 

9 

822.6 

10 

841.7 

11 

880.3 

12 

909.8 

13 

939.8 

14 

969.2 

15 

1000.0 


MEDIUM RANGE 


Meter 

Recorder 

(°C) 

(mv d-c) 

10 

0.0 

11 

24.5 

12 

49.4 

13 

74.9 

14 

100.9 

15 

127.3 

16 

154.3 

17 

181.7 

18 

209.8 

19 

237. 1 

20 

264.8 

21 

292.7 

22 

320.8 

23 

352.0 

24 

377.8 

25 

406.8 

26 

436.0 

27 

465.5 

28 

495.1 

29 

525.0 

30 

554.8 

31 

584.7 

32 

614.2 

33 

644.5 

34 

674.2 

35 

704.8 

36 

735.0 

37 

765.2 

38 

795.7 

39 

829.0 

40 

855.0 

41 

883.5 

42 

912.7 

43 

941.1 

44 

970.1 

45 

1000.0 


HIGH RANGE 


Meter 

Recorder 

(°C) 

(mv d-c) 

40 

0.0 

41 

25. 3 

42 

51.1 

43 

77.3 

44 

104.0 

45 

131.2 

46 

158.8 

47 

186.9 

48 

215.5 

49 

243,4 

50 

271.4 

51 

299.7 

52 

328.2 

53 

356.9 

54 

385.7 

55 

414.8 

56 

441.1 

57 

473.7 

58 

503.5 

59 

535.3 

60 

562.4 

61 

592.1 

62 

622.0 

63 

651.3 

64 

680.8 

65 

711.0 

66 

742.1 

67 

771.1 

68 

801.0 

69 

830.0 

70 

857.6 

71 

885.7 

72 

914.9 

73 

943.0 

74 

971.0 

75 

1000.0 


Note: High-level output is Usted. Low-level output is 5% of value listed for 

each temperature. 


APPENDIX 4 


GRAPHICAL REPRESENTATION OF MEASURED DIURNAL SURFACE 


TEMPERATURES v. THERMAL MODEL RUNS 


SITE 1-1 
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APPENDIX 5 


DETAILED ANALYSIS OF STATION 1-1 EMPLOYING SURTEMP 
THERMAL MODEL - "INPUT TE„p" MEASURED AT SITE 


MARCH 29-30, 1975 
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16128 Trjj 39 HAY 75 


PtOe 3 


<EWBUDI>3|M,|1 


□OOQDOOQOO 

TIMES# FINAL SURFACE TEMPERATURES, AMD INPUT TEMPERATURES 
} 2 »$ 35,1 13,2 27,1 28 » J 27,0 27,7 15,6 24.1 24,5 16,8 18,7 28,1 

l !*S J J #3 19,2 7,5 8,0 20 * a a * 6 3 » 7 21.6 2.4 1,3 22,8 0,7 M.0 

2 ?«J *?•! ’ 2, ° 1,2 " 2 »° * 2 * 3 2 » 4 " 3 »° *2,6 3,6 -4,0 -3,7 4,8 *4^6 *4,4 

6,0 *4,8 -5.4 7,2 -0,8 -2,8 8,4 5,7 4,0 9,6 12,9 U,4 10,8 19.5 10,3 

T MAX 27,1 T MIN -4,8 AVERAGE TEMP 8,32 

PROGRAM PARAMETERS CURRENT VALUES 

2'° C63 ALB 0,365 (7) EMS 0,9850 (8) T IN 0,045 (9) LAT 34,7 
(10) DEC 3,0 (11) DIP 0,0 (12) STR 0,0 (13) TN 230,0 (14) TD 240.0 

PLOT OF FINAL TEMP ( + ) AND INPUT TEMP (*) AGAINST TIME 
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aaaaoaaoDa 

TIMES, FINAL SURFACE TEMPERATURES, AMD INPUT TEMPERATURES 

12.0 22.9 25,1 13,2 25,0 28,1 14,0 24,8 27,7 15,6 22.4 24,5 16,8 18,1 20, * 

18.0 12,6 10,3 19,2 9,6 8,0 20.4 7.6 3,7 21,6 6.0 1,3 22.8 4,8 *1,0 

24.0 3.8 .2.0 1.2 2,9 .2,3 2.4 2.2 -2,6 3.6 1.5 .3,7 4,8 0,9 * 8,8 

6,0 0.4 .5,0 7,2 3,3 .2.8 8.4 8,2 4,0 9.6 13,8 11,4 10,8 19,0 10,3 


T MAX 25,0 T MIN 0.4 AVERAGE TEMP 10,50 


PROGRAM PARAMETERS CURRENT VALUES 

(5) CLOUD 0.2 (6) ALB 0,365 (7) EMS 0.9850 (8) T IN 0.045 (9) LAT 

(10) DEC 3.0 (11) DIP 0,0 (12) STR 0,0 (IS) TN 250,0 (14) TD 


34,7 

250.0 


PLOT OF FINAL TEMP (♦) AND INPUT TEMP (*) AGAINST TIME 
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TIMES, FINAL SURFACE TEMPERATURES, AND INPUT TEMPERATURES 

12.0 30.8 25,1 13.2 32,8 28,1 1-4.4 31,1 27,7 15,6 ?S,7 P4,S 16,8 17,0 20,1 

18.0 6.8 14,3 19.2 5.3 8.0 20,4 3,7 3,7 21,6 2,4 1,3 22,8 1,3 *1,0 

24.0 0,5 *2,0 1,2 *0,3 -2,3 2.4 -0.9 -2,6 3.6 -1,4 *3.7 4,8 -1,J **.4 

6,0 -6,1 -5,4 7,2 -1.4 -2.8 8.4 7.4 4,0 9.6 16,9 11,4 10,8 25.2 20,3 

T MAX 32,8 T MIN -6,1 AVERAGE TEMP 9,75 

PROGRAM PARAMETERS CURRENT VALUES 

(5) CLOUD 0.2 (6) ALB 0,365 (7) EMS 0,9850 (8) T IN 0.020 (9) LAT 34,7 

(10) DEC 3.0 (11) DIP 0,0 (12) STR 0,0 (13) TN 260.0 (14) TD *40,0 


PLOT OF FINAL TEMP (♦) AND INPUT TEMP (*) AGAINST TI M E 
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. 

¥ 

TIMES, FINAL SURFACE TEMPERATURES, AND INPUT TEMPERATURES 

12.0 37.6 25,1 13.2 40,5 28,1 14,4 38.6 27,7 15,6 32.2 24,5 16,8 21,8 *0,1 

18.0 9,4 14,3 19,2 1,8 8,0 20,4 -2,7 3,7 21,6 •5,8 1,3 22,8 -8,2 -1,0 

24,0-10,1 -2,0 1,2-11,7 -2,3 2,4-13.0 -2,6 3,6-14,2 -3,7 4,8-15,2 -4,4 

6,0-14,5 -5,4 7,2 -6,0 -2.8 8,4 6,5 4,0 9,6 19,5 11,4 10,8 30,4 10,1 


T MAX 40 

,5 T MIN -15,2 AVERAGE TEMP 6,84 

\ , 

PROGRAM PARAMETERS CURRENT VALUES 

(5) CLOUD 0,0 (6) ALB 0,365 (7) EMS 0,9850 (8) T IN 0,020 (9) LAT 34,7 

CIO) DEC 3.0 (11) DIP 0,0 C 123 5TR 0,0 C IS) TN 230,0 C14) TD *40,0 


PLOT OF FINAL TEMP (♦) AND INPUT TEMP (*) 
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AGAINST TIME 
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TIMES, FINAL SURFACE TEMPERATURES, 

AND 

INPUT 

TEMPERATURES 
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T MAX 35,6 T MIN -6,9 AVERAGE TEMP 9,69 


PROGRAM PARAMETERS CURRENT VALUES 

(5) CLOUD 0,2 (6) ALB 0,365 (7) EMS 0.965Q (8) T IN 0,020 (9) LAT 34,7 

(10) DEC 3.0 (11) DIP 0,0 (12) STR 0,0 (13) TN 230,0 (14) TD |50,0 


plot OF 

FINAL TEMP (♦) 

AND 

INPUT 

TEMP 

(*) 

AF-AINST TIME 
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TIMES, FINAL SURFACE TEMPERATURES, AND INPUT TEMPERATURES 

12.0 24.8 25.1 13.2 27,0 28,1 14.4 25.4 27,7 15,6 19.9 24.5 16.8 11,1 20. 1 

18.0 0.7 14.3 19.2 -1.3 8,0 20.4 -3,1 3,7 21.6 -4.7 1.3 22.8 -5,9 -1.0 

24.0 -6.9 -2.0 1.2 -7,7 -2,3 2.4 -8.4 -2.6 3.6 -9.0 -3.7 4,8 -9,6 *4.4 

6,0-13,4 -5,4 7,2 -8,6 -2.8 8.4 0,5 4,0 9,6 10,4 11,4 10,8 18.9 20,3 

T MAX 27,0 T MIN -13,4 AVERAGE TEMP 3,01 

PROGRAM PARAMETERS CURRENT VALUES 

(5) CLOUD 0.2 (63 ALB 0,365 (73 EMS 0.9850 (8) T IN 0.020 LAT 14.7 

(10) DEC 3.0 (11) DIP 0,0 (12) STR 0,0 (13) TN 250,0 (14) TD 230,0 


PLOT OF FINAL TEMP (+) AND INPUT TEMP (*) AGAINST TIME 
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TIMES, FINAL SURFACE TEMPERATURES, AND INPUT TEMPERATURES 

12.0 28,9 25,1 13,2 31,1 28,1 14,4 29,5 27,7 15,6 24,2 24,5 16,8 IS*, 5 10,1 

18.0 5.3 14,3 19,2 2,0 8.0 20,4 *0.5 3,7 21,6 *2,3 1,3 22,8 -3.7 -1,0 

24.0 -4,9 -2,0 1.2 -5,8 -2,3 2,4 -6,7 -2,6 3,4 -7,4 -3,7 4,8 -8,0 -4,4 

6,0-10,3 -5,4 7.2 -4,8 -2.8 8,4 4,5 4,0 9,6 14,5 11,4 10,6 23.1 *0,3 

T MAX 31,1 T MIN -10,3 AVERAGE TEMP 6,21 

PROGRAM PARAMETERS CURRENT VALUES 

(5) CLOUO 0,2 (6) ALB 0,365 (7) EM3 0.9850 (8) T IN 0,020 (9) LAT 34,7 

(10) OEC 3.0 (11) DIP 0,0 (12) STR 0.0 (13) TN 250.0 (14) TD 240,0 


PLOT OF FINAL TEMP (♦) AND INPUT TEMP (*> AGAINST TIME 
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TIMES, FINAL SURFACE TEMPERATURES, AND INPUT TEMPERATURES . t 

12.0 37,0 25,1 13.2 39,9 28.1 14.4 38,0 27,7 15,6 31,6 24,5 16,8 21,2 20,1 

18.0 8,8 14,3 19,2 1,3 8,0 20,4 -3,2 3,7 21,6 -6,3 1.3 22,8 -8,7 -1.0 

24.0-10.6 -2.0 1.2-12,2 -2,3 2.4.13,5 -2,6 3.6-14.7 -3,7 4, 8-15,7 -4,4 

6,0-14,9 -5,4 7,2 -6,5 -2,8 8.4 6,1 4,0 9,6 19,0 11.4 10,8 29,8 20,3 

T MAX 39,9 T MIN -15,7 AVERAGE TEMP 6,31 

PROGRAM PARAMETERS CURRENT VALUES . . 

(55 CLOUD 0,0 (6) ALB 0,365 (7) EMS 1,0000 (8) T IN 0.020 (9) LAT 34,7 

(10) DEC 3,0 (11) DIP 0,0 (12) STR 0.0 (13) TN 230.0 (14) TD 240.0 


PLOT OF FINAL TEMP (♦) AND INPUT TEMP (*) AGAINST TIME 
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TAUTOCHRONES FROM PLAYA PROBE MEASUREMENTS 
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